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A HIGHER-ORDER LARGE-SCALE REGULARITY THEORY 
FOR RANDOM ELLIPTIC OPERATORS 

JULIAN FISCHER AND FELIX OTTO 


Abstract. We develop a large-scale regularity theory of higher order for 
divergence-form elliptic equations with heterogeneous coefficient fields a in the 
context of stochastic homogenization. The large-scale regularity of a-harmonic 
functions is encoded by Liouville principles: The space of a-harmonic functions 
that grow at most like a polynomial of degree k has the same dimension as in 
the constant-coefficient case. This result can be seen as the qualitative side 
of a large-scale C k,a -regularity theory, which in the present work is developed 
in the form of a corresponding C k5 “-“excess decay” estimate: For a given a- 
harmonic function u on a ball Br, its energy distance on some ball B r to the 
above space of a-harmonic functions that grow at most like a polynomial of 
degree k has the natural decay in the radius r above some minimal radius r$. 

Though motivated by stochastic homogenization, the contribution of this 
paper is of purely deterministic nature: We work under the assumption that 
for the given realization a of the coefficient field, the couple (0, ct) of scalar and 
vector potentials of the harmonic coordinates, where (f) is the usual corrector, 
grows sublinearly in a mildly quantified way. We then construct “fcth-order 
correctors” and thereby the space of a-harmonic functions that grow at most 
like a polynomial of degree k, establish the above excess decay and then the 
corresponding Liouville principle. 


1. Introduction 

We are interested in the regularity of harmonic functions u associated with a 
uniformly elliptic coefficient field a in d space dimensions (by which we understand 
a tensor field satisfying A|£| 2 < £ • a£ and |a£| < |£| for some A > 0 and any £ £ R d ) 
via the divergence-form equation 

(1) — V- oVw = 0. 

Without continuity assumptions, the local regularity of (weak finite-energy) solu¬ 
tions can be rather low, in particular in case of systems (cf. e.g. [101 Example 3] 
for the scalar case and B Section 9.1.1] for De Giorgi’s celebrated counterexample 
in the systems case). Because of their homogeneity, the same examples show that 
even when the coefficients are uniformly locally smooth, the large-scale behavior 
of a-harmonic functions can be very different from the constant coefficient, that is, 
Euclidean case; cf. e.g. Proposition [5T] in the appendix below. Large-scale regular¬ 
ity is most compactly encoded in a Liouville statement of the following form: The 
space of a-harmonic functions u of growth not larger than \x\ k has the same di¬ 
mension as in the constant-coefficient case, where the space is spanned by spherical 
harmonics up to order k. Because of the above-mentioned counterexamples, such 
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Liouvillo statements may fail for uniformly elliptic coefficient fields: For example, 
in the case of systems, there are non-constant harmonic maps that decay to zero 
at infinity. 

The question whether this situation generically improves for certain ensembles of 
coefficient fields, namely stationary and ergodic ensembles as in stochastic homog¬ 
enization, seems to have first been phrased and partially answered in [4j Chapter 
6 and Theorem 3] in the context of random walks in random environments: Under 
the mere assumption of ergodicity and stationarity, sublinearly growing a-harmonic 
functions are almost surely constant. The argument is limited to the scalar case, 
but can deal with non-uniformly elliptic cases as percolation. 

Motivated by error estimates in stochastic homogenization, the topic of a regu¬ 
larity theory for random elliptic operators was independently addressed in a more 
quantitative way in [9j. In Corollary 4 of that paper, for any a < 1, a large-scale 
C°’ a -inner regularity estimate for a-harmonic functions has been established, with 
a random constant of finite algebraic moments — however under stronger assump¬ 
tions on the ergodicity, namely a finite spectral gap w. r. t. Glauber dynamics in 
the case of a discrete medium. 

A major step forward constitutes [2], where the above result was improved to a 
large-scale C' 0,1 -inner regularity estimate even in case of (symmetric) systems, by 
showing that the approach of [3| for obtaining (large-scale) regularity of a-harmonic 
maps, which itself is based on a Campanato-type iteration, can be extended from 
periodic to random coefficient fields. Under a strong assumption of ergodicity, 
namely that of a finite range of dependence, optimal exponential moments for the 
random constant are obtained. 

This work motivated [7], which in turn is the basis for the present paper. In that 
work, another tool from periodic homogenization, namely the vector potential er for 
the harmonic coordinates (next to the well-known scalar potential (j >, also called the 
corrector), was transferred to the random case, see (0 and 0 for the characterizing 
properties. This allowed to establish a C 1, “-Liouville theorem, meaning that the 
space of sub-quadratically growing a-harmonic functions is almost surely spanned 
by the the constants and the d a-harmonic coordinates Xi+fa. This holds even for 
non-symmetric systems and was shown under the mere assumptions of stationarity 
and ergodicity. More precisely, it relied on the almost-sure sublinear growth of the 
couple (</>, a) of correctors in the sense of 

(2) lim e r = 0, 

r—>oo 

where 

(3) e r := sup -j; (-f \(/)\ 2 + \a\ 2 dx 

R>r H \Jb r 

This sublinear growth 0 was shown to hold under the assumptions of stationarity 
and qualitative ergodicity. In a second step, large-scale C 1 ’“-inner regularity esti¬ 
mates for a-harmonic functions were obtained, where the random constant satisfies 
a stretched exponential bound under mild decay assumptions on the spatial covari¬ 
ance of a. In a later version of [7] , the optimal stochastic moments for the random 
constant were obtained. 

In the context of non-linear elliptic systems in divergence form, the result of [2] 
on the large-scale C' 0,1 -estimate was generalized in pQ to non-symmetric coefficients 
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and well beyond finite range, further confirming that the random large-scale regu¬ 
larity theory holds under just a mild quantification of ergodicity, like expressed by 
standard mixing conditions. 

In the present work, we go beyond C 1,a and establish a large-scale C k,a -theory 
in form of a corresponding excess decay and Liouville result, cf. Theorem [3] and 
Corollary [4] This lifts the result of [S’ from the periodic to the random case. To 
streamline presentation, we first establish the C 2 ’“-versions of our theorems, cf. 
Theorem [7] and Corollary [5] 

Let us clearly state that the contribution of this paper is exclusively on the 
deterministic side. The large-scale regularity is obtained under the assumption 
that the given realization a of the coefficient field is such that the corresponding 
corrector couple (<p, a) satisfies the following slight quantification of ([2]), namely 

(4) lim S2r = 0 
with 

OO 

(5) e 2 , r ■= ^2 m ^ n {l, 2 m+1 /r}e 2 m. 

m=0 


Note that © is equivalent to X^m=o e 2 m < 00 ■ 

In the recent preprint [5], it is shown that © holds for almost every realization 
a in case of a stationary ensemble of coefficient fields under mild quantification of 
ergodicity in form of an assumption on a mild decay of correlations of a: More 
precisely, given a stationary centered tensor-valued Gaussian random field a on 
R d and a bounded Lipschitz map $ : R dxd —> R dxd taking values in the set of 
A-uniformly elliptic tensors, the coefficient field a := 4>(d) almost surely admits 
correctors with the property © assuming just decay of correlations in the sense 
|(a(:r)d(y)}| < C\x — y | -/3 for some C > 0 and some /3 G (0,c(d, A)) (where (•) 
denotes the expectation). Note that under the assumption of a spectral gap for the 
ensemble, as far as the corrector cp is concerned (but not the “vector potential” cr), 
an estimate like © could also be deduced to hold almost surely from j5j Proposition 
2], modulo the passage from a discrete to a continuum medium. 

The key building block for this large-scale C fe, “-theory is the space of a-harmonic 
functions that grow at most like a polynomial of degree k at infinity. Proposition 
[2] and Corollary |4] imply that under our assumption © this space has the same 
dimension as in the Euclidean case - e.g. for k = 2 the space of a-harmonic functions 
that grow at most quadratically is spanned by l+d+ d ^ ^ maps - , which partially 
answers the question in TJ Chapter 6]. The fcth-order excess (fill) , by the decay 
of which we encode the C fc,a -theory, measures the distance to this space in terms 
of the averaged squared gradient. As our construction shows, there is a one-to- 
one correspondence between the asymptotic behavior of functions in this space 
and aft, om -harmonic polynomials of degree k. However, there is no natural one-to- 
one correspondence between elements of this space and /cth-order a^om-harmonic 
polynomials. 

Before stating our results, let us recall the definition of the correctors ( cp , a). The 
corrector (pi satisfies the equation 

(6) - V-a(ei + Vcpi) = 0. 
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The flux correction qij is defined as 

(7) qi • = d{Ci T dhomE-i 

where dhom is the homogenized tensor, that is, dhom^i is the expectation of a(ej + 
V0i). In our analysis, we will only use that dhom is some constant elliptic coefficient. 
We introduce the corresponding vector potential a l3 k (antisymmetric in its last two 
indices) by requiring that 

(8) V- <Jij = qij. 

For the actual construction of a a with stationary gradient we refer to [7]; in this 
note, we just use the property ((SI). In the context of periodic homogenization, 
both the scalar and the vector potential <p and a may be chosen to be periodic. 
In stochastic homogenization, one cannot always expect to have a stationary (if), a) 
(for instance in d < 2 even in case of finite range) but, as mentioned above, we 
expect sublinear growth in the sense of ([4]) under mild ergodicity assumptions. 

Notation. Throughout the paper, we use the Einstein summation convention, 
i.e. we implicitly take the sum over an index whenever this index occurs twice. For 
example, bidiV is an alternative notation for (b ■ V)v and biS7vi is an alternative 
notation for Xu=i biVvi. 

By C we denote a generic constant whose value may be different in each appear¬ 
ance of the expression C; similarly, by e.g. C(d, A) we denote a generic constant 
depending only on d and A whose value again may be different for every use of the 
expression C(d, A). 

By £ := {E £ R dxd : (Eij + Eji)(dhom)ij = 0} we denote the space of matrices 
Eij for which EijXiXj is an dhom -harmonic second-order polynomial. 

The notation P (or P(x)) generally refers to a polynomial. By V k , we denote the 
space of homogeneous polynomials of degree k. By V k hom , we denote the space of 
homogeneous polynomials of degree k which are afc om -harmonic. On the space V k , 
we introduce the norm ||P|| := sup a . gBl |P(a;)|; note that any other norm on this 
finite-dimensional space would do as well, since we do not care for C(fc)-constants. 

2. Main Results 

The proof of our large-scale C k,a regularity theory relies in an essential way on 
the existence of fcth-order correctors for the homogenization problem, which enable 
us to correct a^om-harmonic polynomials of degree k by adding a small (in the 
L 2 -sense) perturbation. 

The ansatz for the deformation of an a^om-harrnonic polynomial P, homogeneous 
of degree k (i.e. P E 'Pa hom )> into an a-harmonic function u with the same growth 
behavior is motivated by homogenization: We consider P as the “homogenized 
solution of the problem solved by u”, so that we think in terms of the two-scale 
expansion u « P + (j>kdkP and have that the error ipp := u — (P + (f>kdkP ) satisfies 
—V • aVf/ip = V • ((cf>kd — Ofc)V<9fcP). In order to construct u, we reverse the 
logic and first construct a solution ipp to the above elliptic equation and then set 
u := P + (f> k d k P + ipp. 

Theorem 1 (Existence of higher-order “correctors for polynomials”). Let d > 2, 
k > 2, and suppose that the corrector <p and the flux-correction potential a satisfy the 
growth assumption 0- Letro be large enough so that e- 2 , ro < £o holds (the existence 
of such ro is ensured by m where Eq = Eo(d,k,X) > 0 is a constant defined in 
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the proof below. Given any P G V k , there exists a “corrector for polynomials’’ ipp 
satisfying 

(9) — V-aVi/jp = V-((^ja — 

as well as 

(10) ^p_L_(£ \V,p P \ 2 dx)j <C(d ) fc 1 A)||P||e 2 , r 
for any r > tq- Moreover, ipp depends linearly on P. 


Our xpp indeed enable us - in conjunction with the first-order correctors (pi - to 
correct a/j 0m -harmonic fcth-order polynomials. 

Proposition 2. Let d > 2, k > 2, and let P G P* . Suppose that ipp satisfies 
©• We then have 

— V- aV (P + (pidiP + ipp) = 0. 


Let us now state our C k,a large-scale regularity result. 

Theorem 3 (C fc ’“ large-scale excess-decay estimate). Let d > 2, k > 2, and sup¬ 
pose that © holds. Let u be an a-harmonic function. Let ipp = 0 for linear poly¬ 
nomials P (in order to simplify notation) and let ipp be the functions constructed 
in Theorem\]\for higher-order polynomials. Consider the kth-order excess 

2 


(11) ExCfc(r) := inf 


p K ev?. 


k 


Vu - V ^(P K + (pidiP K + ipp K ) 


K = 1 


dx. 


Let 0 < a < 1 and let r$ be large enough so that £2 , ro < £0 holds (the existence of 
such ro is ensured by ^)), where £0 = £ 0 (d, k, A, a) > 0 is a constant defined in the 
proof below. Then for all r, R > ro with r < R the C k,a excess-decay estimate 

/ T \ 2(fc—l)+2ai 

Exc fc (r) < C(d, k, A, a) (— ) 


(12) Exc fc (r) < C(d, k, A, a) J Exc k (R) 

is satisfied. 

Our large-scale C k+1,a excess-decay estimate entails the following fcth-order Li- 
ouville principle. 


Corollary 4 (fcth-order Liouville principle). Let d > 2, k > 2, and suppose that 
the assumption © is satisfied. Then the following property holds: Any a-harmonic 
function u satisfying the growth condition 

1/2 

(13) 

is of the form 


liminf ^ ( / |u| 2 dx ) =0 

r—f 00 y'k-rL \ l g 


u = a + bi(xi + (pi) + ^(P K + (pidiP K + i/tpj 


K =2 


with some a€l, t 6 R d , and P K G "Pa hom f or 2 < k < k (i.e. P K is a homogeneous 
CLhom-harmonic polynomial of degree k). Here, the ipp denote the higher-order cor¬ 
rectors whose existence is guaranteed by Theorem. QJ 

In particular, the space of all a-harmonic functions satisfying da has the same 
dimension as if a was replaced by a constant coefficient, say ahom- 
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The structure of our proofs is as follows: 

• First, under the assumption that we already have constructed an appropri¬ 
ate fcth-order corrector on a ball Br, we show a C k,a excess-decay estimate 
on large scales within this ball for a-harmonic functions (Lemma H4l) . This 
result directly implies Theorem [3] as soon as we have proven the existence 
of a corrector on W 1 (i.e. Theorem [lj. 

• Our C k,a estimate implies a theory for a-harmonic functions on balls 

Br, provided that we have already constructed an appropriate fcth-order 
corrector on Br. This is done in Lemma [171 

• At last, we are able to build our corrector, starting from small balls and 
iteratively doubling the size of our balls. At this point, we require the 
C fc_1 ’ 1 theory to show appropriate (fcth-order) decay in the interior of the 
new contribution to the fcth-order corrector entering at every scale. This 
iterative enlargement is carried out in Lemma [18] and finally enables us to 
prove Theorem [T] 

• The fcth-order Liouville principle stated in Corollary [I] is an easy conse¬ 
quence of our C' fc+1, “ large-scale excess-decay estimate. 

3. A C 2,a Large-Scale Regularity Theory for Homogeneous Elliptic 
Equations with Random Coefficients 

For the reader’s convenience, we shall first provide a proof for the C 2,a case of 
our theorems, as in this case the proofs are less technical while already containing 
the key ideas. In particular, the overall structure of our proofs is the same as in 
the C k,a case. Since we shall use a somewhat simplified notation in the C 2,a case, 
let us reformulate the C 2,a case of our theorems using this notation. 

Theorem 5 (Existence of second-order correctors). Let d > 2 and suppose that the 
corrector <p and the flux-correction potential cr satisfy the growth assumption 0- 
Let ro be large enough so that £ 2 ,r 0 < £o holds (the existence of such vq is ensured 
h where Eo = Eo{d , A) > 0 is a constant defined in the proof below. Given any 
E £ R dxd , there exists a second-order corrector ip E satisfying 

(14) - V- aS/ipE = Eij V-[<r y + crji + afaej + <^e»)] 

as well as 

su P i(7 \Vtp E \ 2 dx) <C{d 1 X)\E\£ 2 ,r 

R>r E \Jb r / 

for any r > ro- Moreover, the corrector S7ip E depends linearly on E. 

Due to the linear dependence of ipE on E, below we shall also write Eijipij in 
place of ipE- 

Note that our second-order correctors indeed enable us - in conjunction with the 
first-order correctors (pi - to correct a^ OTO -harmonic second-order polynomials. 

Proposition 6. Let d > 2 and let E £ £ (i.e. assume that the polynomial EijXiXj 
is ahom-harmonic). Suppose that ip e satisfies (fl4l) . We then have 

— V- aVEij(xiXj + Xi<pj + <piXj + tpij) = 0. 

Our C 2, “ large-scale regularity theorem reads as follows. 
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Theorem 7 (C 2,a large-scale excess-decay estimate). Let d > 2 and suppose that 
© holds. Let u be an a-harmonic function. Let i[)e be the second-order corrector 
constructed in Theorem [3 Consider the second-order excess 

(15) 

Exc 2 (r) := 

inf / \\/u-\/(bi(xi + (f>i) + EiAxiXj + Xi(/)j + (fiXj+4> ii ))\ 2 dx. 

bm d ,EG£j Br 1 v n 

Let 0 < a < 1 and let ro be large enough so that e 2 , ro < £o holds (t/ie existence of 
such ro is ensured by m where Eq = £o{d, A, a) > 0 is a constant defined in the 
proof below. Then for all r 1 R > ro with r < R the C 2,a excess-decay estimate 

(16) Exc 2 (r) < C(d, A, a) J Exc 2 (E) 
is satisfied. 


Our large-scale excess-decay estimate entails the following C 2,a 
pie. 


Liouville princi- 


Corollary 8 ( C 2 ’ a Liouville principle). Let d > 2 and suppose that the assumption 
© is satisfied. Then the following property holds: Any a-harmonic function u 
satisfying the growth condition 


lim inf 


1 

~\~ot 




= 0 


for some a £ (0,1) is of the form 


u = a + bi(xi + (fi) + Eij(xiXj + Xirfij + cfiXj + ifij) 


with some a £ R, b £ and E £ £ (i.e. some E £ WL dxd for which EijXiXj is an 
ahom-harmonic polynomial). 


Let us start with the proof of Proposition [6j which only requires a simple com¬ 
putation. 


Proof of Proposition [6J Making use of the fact that Eij((ahom)ij + ( Phom)ji ) = 0 
(in the third step below), we compute 

Ejj V-(cq^ T c 7ji ) T Eij V' a((f e.j T cfjefj 

— dij ~L Ej j Qji Eij V - O ( 0i Cj T (f)j Cj ) 

= Eij (ttjTc ( (1 d ) I/c “t“ ( ahom)ji ) T Elj (Aik ((Id)y'A: T 9k<pj) (ahom'jij) 

E/j V* of. 0iCj ~\- 0jCl ) 

— Epy (Ujfc (fikXi T Of .Oj ) T UzA: {Of-Xj T Of.(f)j )) 

+ Eij V- a((f>iVxj + (fjS7xi) 

= Eij(aV(xi + (j>i) ■ S7xj + aV(xj + </>j) ■ Vxi) + E t j V- a(c/) t X7xj + ^Vs,) 

= Ey V-(xjaX7(xi + <pi) + XiaV(xj + (f>j)) + E^ V- a((f>iX7xj + (f>jS7xi). 
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We therefore obtain 


Eij V-(cTy + (Tji) + E i:j V- a(4>iej + erfj) 

= V- aV(xiXj + Xi4>j + fcxj), 

which together with (1141) implies our proposition. □ 


3.1. The C 2, “ excess-decay estimate. To establish our C 2,a excess-decay esti¬ 
mate, we make use of the following lemma, which essentially generalizes Theorem [7] 
to correctors which are only available on balls Br. 


Lemma 9. Let d > 2. For any E £ £, denote by ipE a solution to the equation 
of the second-order corrector (ED on the ball Br (without boundary conditions); 
assume that i/jr depends linearly on E. Set 

( 17 ) £ ti,r,R : = SU P P" 1 max / \V^ E \ 2 dx 

r<p<R y.E6£,|.E| — lj R p 

For an a-harmonic function u in Br, consider the second-order excess 



(18) Exc2(r) := 


inf / 

beR d ,E££J B r 


Vu - V (bi{xi + 4>i ) + Eij[xiXj + Xi<t>j + 4>iXj + -i lij)) 


dx. 


For any 0 < a < 1 there exists a constant £ m in > 0 depending only on d, X, and 
a such that the following assertion holds: 

Suppose that ro > 0 satisfies e ro + ro r < e min- Then for all r £ [ro, R] the 
C 2,a excess-decay estimate 

-—— / v \ 2 + 2 o: - 

(19) Exc 2 (r) < C(d,X,a) Exc 2 (i?) 

is satisfied. 

Note that the infimum in ra is actually attained, as the average integral in 
the definition o/Exc 2 (p) is a quadratic functional of b and E. Denote by 5 p,mm 
and E p,mm a corresponding optimal choice of b and E in We then have the 

estimates 


( 20 ) 


R 2 \E r ’ min - E R ' min | 2 + | b r ’ min - b R ’ min | 2 < C(d, A, a)Exc 2 (i?) 


and 

(21) i? 2 |E I '’ min | 2 +|6 r ’ min | 2 < C(d,A,a)/ \Vu\ 2 dx. 

Jb r 


Proof of Theorem [?] Theorem [7] obviously follows from Lemma [9] by setting %pE '■= 
ipE, with i(e being the second-order corrector whose existence is guaranteed by 
Theorem [5j □ 


The following lemma is essentially a special case of our C 2,a large-scale excess- 
decay estimate Lemma O it entails the general case of Lemma [9] cf. below. 

Lemma 10. Let d > 2 and let R, r > 0 satisfy r < R /4 and £r < 1. For any 
E £ £, denote by ipE a solution to the equation of the second-order corrector 
on the ball Br (without boundary conditions); assume that ipE depends linearly on 
E. For an a-harmonic function u in Br, consider again the second-order excess 
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m- Then the excess on the smaller ball B r is estimated in terms of the excess on 
the larger ball Bn and our quantities £r and Vf/'B' We have 


Exc2(r) < C(d, A) 


(£)' 


2/(d+l) 2 
+ 1 t R 


3-2 


It ~ max 

E££,\E\ = iJ Br 




x Exc 2 (-R). 

Before proving Lemma fTUl we would like to show how it implies Lemma [9] 

Proof of Lemma [3 First choose 0 < 6 < 1/4 so small that the strict inequality 
C(d,X)9 4 < Q 2 + 2a is satisfied (with C(d, A) being the constant from Lemma fTUl) . 


Then, choose the threshold £■ 


mm for £ ro + £^ rg R so small that the estimate 


C(d, A) 


9 4 + I e. 


2 /(d+l) ; 

ro 


+ £ 2 j A 6 - 

■0,r o ,-R/ 


^ q2-\-2c 


holds. 

Let M be the largest integer for which 0 M R > r holds. Applying Lemma [TUI 
inductively with R m := 9 m ~ 1 R 1 r m := 9 m R for 1 < to < M, we infer 

Exc 2 (0 M i?) < (9 2+2a ) M F^e 2 {R). 

Since we have trivially 


Exc 2 (r) < ( — ) Exc 2 (r M ) 

\r M ; 


—d 


and since by definition of M we have r > Qvm and thus 0 M <6 (where we 
recall 9 = 9(d , A, a )), we infer 

Exc2(r) < C(d,\,a ) ^ — J Exc2(i?). 

It remains to show the estimates for | b r ’ mm — b R ' mm \ and | E r ' min — E R ' mm \ as well 
as the bounds for |& r ' ,rram | and \E r,mm \. To do so, let us first estimate the differences 

_ b r m ,rnin | &nd | E R m ,min _ Wg havg the est i mat e 


L 


IV (fo. 


R m , min 7 r m , min 


X )(Xi + <fi) 




Rm,min 
ij ~ ^ij 


)(xiXj + Xifij + cfiXj + I dx 


< 2 / 

Jb, 


IB, 
+ 2 


Vu - V+ 


(xi + <j>i) - V Elp’ mm (xiXj + Xi<j>j + tj>iXj + tyij) 


dx 


t 


Vu — V6 


Rm ,min 


< 2Exc 2 (r m ) + 2 


Rn 


(x-i + cf>i) -VE t 


LXC 2 ( Rm ) 


Rm ,min 

ij 


(xiXj + Xi4>j + 4>iXj + Ipij) 


dx 


< 


< 


, r \2+2a _ , (R \ z + za _ 

C{d, A, a) (-g) Exc 2 (i?) + C(d, A, a)9~ d ( -jj- J Exc 2 (I?) 
C(d,X,a)^Y (9 2a ) m: E^c 2 (R). 


From Lemma ITT1 below, we thus obtain 


I b 


R m ,min _ yr m ,min 


+ R\E 


Rm ,min 


-E 


1 Tm ,rmn 1 


< C(d,A,a)(0 Q ) m V Exc 2+O- 
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Note that a similar estimate for the last increment | for m, mm _ yr,mm | _ 

E r ’ min | can be derived analogously. Taking the sum with respect to m and recalling 
that Ri = R and r m = R m +i, we finally deduce 

M r—— - 

| b R ’ min - b r ' min | + R\E R ' min - E r ’ mm | < C(d, A, a) ^ (d a ) m \j Exc 2 (i?) 

m=0 


< C(d, A, a)y Exc2(i?). 

It only remains to establish the last estimate for |6 r,rram | and \E r,mm \. By the 
previous estimate, it is sufficient to prove the corresponding bound for b Ii ' rnm and 
This in turn is a consequence of the inequality 


£ 


Br 

< 


Vb 


R,min / 


(£j + 4>i) + V E R -’ mm (xiXj + Xi4>j + + ipij) 


dx 


2Exc 2 (R) + 2 f |Vu| 2 dx < 4/ |Vu| 2 dx 

J Br JBr 


together with Lemma [TT] below. 

The following lemma quantifies the linear independence of the corrected polyno¬ 
mials Xi + 4>i , Eij(xiXj + Xi<f>j + tpiXj + it is needed in the previous proof. 

Lemma 11. Suppose that for every E £ £ \ {0}, the functions <f> and ipE satisfy 


□ 


2 / dx + p- 2 \E\~ 2 l IViM 
JB„ J B n 


2 dx < £q, 


where £o = £o(d) is to be defined in the proof below. Then for any b £ R d and any 
E £ £, we have the estimate 

\b\ 2 +p 2 \E \ 2 

(22) < C(d)f |V&j(xj + + VEij(xiXj + Xi<pj + fcxj + ipij)\ 2 dx. 

J B p 

Proof. Poincare’s inequality (with zero mean) and the triangle inequality imply 
(j- \Vbi(xi + fa) + VEij(xiXj + Xifa + <f>iXj + 4>ij )| 2 da:j 


> 


> 


1 1 


inf f | bi{xi + 4>i) + Eij(xiXj + Xi<pj + faxj + %pij) - a| 2 dx 
C{d) p aeR \J Bp ) 

1 ^--7 |6 


\ 1/2 


C{d)p 


1/2 


inf 

aeR 


uiXi + EijXiXj — a| 2 dx 


- inf 

a£R 




1/2 


(pi Ej.j(X{(j)j T (j){Xj T 'ipij') a| dx 


On the one hand, by transversality of constant, linear, and quadratic functions we 
have 

\ 1/2 


- inf 
p a£R \Jg 


f I biXi 

J B n 


+ EijXiXj — a| 2 dx 


> 


1 


C(d) 


p\E\). 
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On the other hand, we have by the triangle inequality and Poincare’s inequality 
/ \ 1/2 

I b t </),; -f- E^j (xj(bj T (f>iXj -f" ^Pij ) 1 

J B„ 


- inf 
p \J]3 


a| 2 dx 


< C(d) 


/ \ 1/2 / \ 1/2 ' 
(\b\ + p\E\)~ l-f \(j)\ 2 dx] + p\E\— max (/ IVi/’gl 2 dx) 

P \Jb p J P ee£,\e\=i \J Bp J 


Putting these estimates together, by boundedness of the integrals in the previous 
line by e^p 2 our assertion is established. □ 

Proof of Lemma 1 10\ In the proof of the lemma, we may assume that 

(23) Ekc 2 (i?) = / |Vu| 2 dx. 

Jb r 

To see this, recall that the infimum in the definition of Exc 2 (i?) is actually attained. 
Denote the corresponding choices of b and E by b mm and E min . Replacing u by 
u — b r f lm {xi + (fi) — E™ m (xiXj + Xif)j + cfiXj + ifij), we see that we may indeed 
assume (EDI) : The new function is also a-harmonic due to © and Proposition [G] 
We then apply Lemma [2D] below to our function u. This yields an ahom~ harmonic 
function Uhom close to u which in particular satisfies 

j- \Vu h om \ 2 dx < C(d,X)S |Vu| 2 dx. 

J B R J 2 J Br 

By inner regularity theory for elliptic equations with constant coefficients, the ahom- 
harmonic function Uhom satisfies 

|Vu fco m(0)| + R SUp \V 2 U ho m\ + R 2 SUp | V*U ho m\ 


B 


R /4 


B 


R/4 


< 


( \ 1/2 

/ \Vu hom \ 2 dx) <C{d,X)(-f |Vit| 2 dx'] . 
Jb r/2 J \Jb r J 


Let us define 

b R,Taylor : = Vu/lorn (0), 

E R,Taylor , = V* Uhom (Q). 

Since — V- ahom^Uhom = 0 holds, we infer E R ’ Taylor (ahom)ij = 0 and therefore 

E R,Taylor g £ ( nQte tha t E ^ Ta V 1 ^ = E R,Taylory g y Taylor , g expansion of 

around i = 0we deduce for any x £ Er/\ the bound 


Vu hom (x) - b R ’ Tayl ° r - E W«yior {xj e i+Xi e j ) 


Making use of the identity 


< \x\ 2 sup |V 3 u? lom |. 


5 R/4 


(Id +(v^) 4 ) Vuftom - V ^b R ’ TaVlm (xi + (fi) + -E R ’ Tayl ° r ( XiXj + Xi<f)j + 4>iXj) 


1 E K,Taylor {(j)jei + (j)iej) 


= (Id +(y<f>) t ) (vu hom (x) - b R ' Taylor - , 
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the previous estimate yields in connection with the bound for |V 3 u/ lom | and r < R/ 4 


L 


(Id+(V</>) t )V?p lom - V ( bf ’ Taylor (xi + (/>i) + -E^’ Taylor (xiXj + Xifa + fax fa 


I E R ’ T °- yl °fafae i + faefa 


2 13 


dx 


<C{d,\)(^] 4 1 \S7u\ 2 dxxl |Id+(V</>) t | 2 dx. 

JBr JB r 

By the Caccioppoli inequality for the a-harmonic function Xi + fa (cf. ©), we have 
(24) J- |Id + (V^)*| 2 dx < f \ X J r( j ) ] 2 dx<C{d,\){l+£\ r ). 

JB r r JB 2 r 

The approximation property of Uhom + fadiUhom in -Br /2 from Lemma [20] below 
implies 


/ |Vu ^('Uhom E (friOiUhom) \ dx^C(d : \)cR 

JBr 

Combining the last three estimates and the equality 


2/(d+l) ; 


(i) 7 j v “i ,<b - 


Vu-V b; 


R,Taylor / 


(Xi + fa) + lE^’ Taylor (xiXj + Xifa + (f>iXj + faj) 


(Id +(V^)‘)VMfc om - V ( bf' Tavlor (xi + fa) + ^ E*’ lavlor (xiXj + Xifa + fax fa 


-iR,Taylor / 


l 


+ Z E «’ Tayl ° fafae i + faefa 


13 


1 


- ^ E % TaVl ° r {<t>3 e i + <^ e i + V ^b) 


+ 


V (Uhom H - tyi^i'U'hom) 

we infer 

Vu - v (. b^ Taylor ( Xi + fa ) + - v ETa y lor 

R,Taylor 


k 


B, 

< 4 


+ dilLfiojTil 


[ L 


(xi + fa) + -EfaJ ay or (xiXj + Xi(j)j + faxj + fafa 


(Id +fa7fa t )Vu hom - V ( bfa laywr (xi + fa) + lE R ’ Taylor ( Xi Xj + Xifa + fax fa 


dx 


Jb, 


+ \ E R fa Tayl ofafae i + faefa 


iR,Taylor ( 

ij 


dx 


^^Pij ) 

+ 4 ~r [Vri V \Uhom 4“ fii&i^hom) \ dx 
JB r 

+ 4L \faS7diU h om\ 2 dx 
J B r 

<C'(d,A)(-^) 4 (l+ £ 2 )£ |Vu| 2 dx 


dx 


'Br 

+ C{d)\E R ’ Taylor \ 2 f r 2 e 2 + max 


EES ,\E\=lJ g 


/ NM 2 

1 J Br 


dx 
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+ C(d, A)4 /(d+1)2 (L) |Vu | 2 dx 

+ C(d)r 2 e 2 sup \V 2 u hom \ 2 . 

Br/4. 


This finally yields in connection with the above bounds on V 2 u? lom in Br/ 4 (recall 
that E R ’ Ta y lor = S7 2 u hom ( 0)) 


i 


V it - V [b 


R,Taylor 


(xi + 4>i) + —E, 


R, Taylor 


(xiXj + XiCpj + 4>iXj + ^ ij) 


C'(d,A)(L) 4 (l +£ 2)£ |Vit | 2 dx 


2 2 , 
r e r + max 


+ C(d, X)R~ 2 -f |Vu | 2 dx 
Jb r 

+ C{d,X)£ 2 jl {d+1) (L) y |Vu | 2 da; 
+ C(d, A)r 2 e 2 i? -2 / |Vu | 2 da; 

J Br 

< C(d, A) [(1) 

x / |Vu| 

J Bn 


E££ ,\E\ = lJ g 


/ |VV’ B | S 
ids,. 


dx 


+ 

2 day 


b r 

2 /(d+l ) 2 


+ i? max 


EG£,|i3|=l7B R 


/ivfep *)(!)■■' 


da; 


where in the last step we have used the inequality e 2 < (j-) d £r < (-y) d ■ 

The new bound directly implies the desired estimate. □ 


3.2. The C 1 ’ 1 excess-decay estimate. We now show how our C 2,a excess-decay 
estimate for the second-order excess Exc 2 from Lemma[9]entails a C 1,1 excess-decay 
estimate for the first-order excess Exc. 


Lemma 12. Let d > 2 and R > 0. For any E £ £, denote by ipE a. solution to 
the equation of the second-order corrector m on the ball Br (without boundary 
conditions); assume that %Pe depends linearly on E. There exists a constant £ m i n > 
0 depending only on d and A such that the following assertion holds: 

Suppose ro £ (0, R] is so large that £ ro < £ m in and 


sup p 

r 0 <p<R 


max 
Ee£,\E\=lJ B 


/ Iv<M 2 

JB n 


1/2 


dx 


< £. 


_ °min 


hold. Let u be an a-harmonic function on Br. Then there exists b R £ for which 
the estimate 



X7b R (xi + (pi)\ 2 dx<C{d,X)(-\ / |Vu | 2 dx 

KRJ Jbr 


holds for any r £ [ro,i?]. Furthermore, b R depends linearly on u and satisfies 


\b R \ 2 < C(d, X)-f \Vu\ 2 dx. 

J Br 
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Proof. In Lemma 0 fix a := 1/2. We then easily verify that Lemma 0 is applicable 
in our situation. Set b R := b r0,mm and E R := E r °’ min ; this implies that b R depends 
linearly on u. The estimate m takes the form 

R 2 \E R \ 2 + \b R \ 2 <C{d,X)-f |Vu | 2 dx. 

Jb r 

Furthermore, applying Lemma0with r 0 playing the role of r and r playing the role 
of R , we deduce from (0(71) 


r 2 \E R - E r ’ min | 2 + | b R - b r ’ min | 2 < C(d, A)Exc 2 (r) 

p \ 2 -|- 2 q ' 


(Tol l 

< 


( p \ z-\-’Zen / p \ 2~\~2oi C 

-J Exc 2 (/i)<C(d,A)(-J |Vu | 2 dx. 


We now estimate 


f \Vu-X/b R {xi+ (j>i)\ 2 dx 
J B r 

l L 


IB, 

< 3 


Vu - Vb r i ' mm (x i + <j>i) - V Ef: j min (xiXj + Xi<j)j + <f)iXj + i>ij) 


dx 


V U[’ mm (x^j + Xi<j)j + <t>iXj + ipij) 


+ 3 / 

JB n 

+ 3/ | (bi’ m%n - 6f )V(®i + fa) | 2 dx 

J B r 


dx 


< 3Exc2(r) 

+ C(d)\E r ’ min \ 2 


/ \(j)\ 2 + r 2 | Id+(V<W | 2 dx + max 
Jb' Ee£,\E\=lJ B 


/ |Vi/e | 2 

Br 


dx 


! (e 2 + (1 + e|r) + £ l,r 0 ,R> 


+ 3| b r ' min - b R \ 2 -f I Id PiVcfYl 2 dx 

J Br 

1191241 / r \ 2+2a__ 

< C(d, A)(-J Exc 2 {R) + C(d 1 X)\E 

+ C(d, X)\b r,min — b R \ 2 (l + s 2 r ) 

< C(d, A) (-J Exc 2 (i?) + C(d, \)\E r ’ min \ 2 r 2 + C(d, \)\b r ’ min - b R \ 2 . 

In conjunction with the two previous estimates, we infer 
/ |Vu - V&f (xi + f>i)\ 2 dx 

J Br 


< C(d, A) 


/ p \ 2+2a ( f P \ 2 / r \ 2+2o;\ / p \ 2+2o:l C 

(b) + ((h) Hr) Mb) \L' Vu ' dx ■ 


Our lemma is therefore established. 


□ 


3.3. Construction of second-order correctors. Using the C 1,1 theory estab¬ 
lished in the previous subsection, we now proceed to the construction of our second- 
order corrector. The following lemma provides the inductive step; starting from a 
function which acts as a corrector on a ball B Bl we construct a function acting as 
a corrector on the ball B 2 r. 
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Lemma 13. Let d > 2 and let ro > 0 satisfy the estimate £2 ,r 0 < £ 0 , where 
£0 = £o{d, A) is to be chosen in the proof below. Then the following implication 
holds: 

Let R = 2 M ro for some M £ No- Suppose that for every E £ R dxd we have a 
solution ip R to the equation 

- V- aVf if = Eij V- XB R [&ij + o-ji + a{cjHej + (pjef)] 


subject to the growth condition 

v. 1/2 M 

|VV’b| 2 dx j < Ci(d,X)\E\ ^2 min{l,2 m r 0 /r}£2“r o 

' m =0 

for all r > ro, where C\{d, A) is a sufficiently large constant to be chosen in the 
proof below. Assume furthermore that ip R depends linearly on E. 

Then for every E £ M. dxd there exists a solution ip^ to the equation 

S7~ a^ip R T cr.j ,■ -f- ofcbjCj 4- cpj Ci ))] 



subject to the growth condition 


/ n \ 1/2 M-\- 1 

r_1 (/ dx ) - C iK'MI- E l X! mi n{ 1 >2 mr oA}£2 

V/B,. / m=0 

/or all r > ro- Furthermore, ip^P depends linearly on E and we have 

' r \ 1/2 

J~ b \Vlp 2 E R - V0 || 2 otej < C^d, A)|E|£ 2 M+i ro . 


Proof. To establish the lemma, we first note that the assumptions of the lemma 
ensure that the C 1,1 excess-decay lemma (Lemma fl2l) is applicable on Bn with 
ipE ■= iPe- To see this, we estimate for any r £ [ro, 1?] 



< C\(d , A)|£f|£2,ro 


< Ci(d, A)j£7|£o- 


By choosing £ 0 > 0 small enough depending only on d and A and C\ (which is to be 
chosen at the end of this proof), we can ensure that the assumption of Lemma l~H/l 
regarding smallness of £7 ro R is satisfied. 

Let now £j[ be the weak solution on R d with square integrable gradient, which is 
unique up to additive constants and whose existence follows from the Lax-Milgram 
theorem, to the problem 

— V- uV/jj = Eij V- XB 2 R—B R (o'ij + 0J») T Eij V- XB 2 R — BR < l{ ( Pi£j + <Pj£i)- 


Obviously, V£§ depends linearly on E; after fixing the additive constant e.g. by 
requiring J B dx = 0, itself depends linearly on E. Furthermore, we have the 
bound 

[ |V£f | 2 dx < C(A)|.E | 2 [ Xb 2 r-BrW\ 2 +Xb 2 r-Br\(P\ 2 dx 
J R d 

and therefore 

(25) [ |V /|| 2 dx < C(\)\E\ 2 R 2+d £ 2 R . 

J R d 
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As£f is a-harmonic in B R , Lemma fl2l now implies the existence of some b E £ 
for which the estimates 


(26) \b%\ 2 < C(d, X)l |V£f | 2 dx < C(d , X)\E\ 2 R 2 s 2 2R 

Jb r 

and 


£ |V££-V(&£)i(*i + &)| 2 dx<C(d,\)(£f f- \V£§\ 2 dx 

<C(d,X)\E\ 2 r 2 e 2 R 


hold for all r £ [ro,-R] and which linearly depends on E. 
Furthermore, we have for r > R 


L 


— S7{bfi)i(xi + (f>i)\ 2 dx 


l24t 

< C(d,A)(r 
1251261 


d 


IB, 

12 D 2 


iveii 2 dx + |6|| 2 (i+ e i) 

d \ 


R 


< C(d,X)\E\ 2 R 2 \ - +1 + 4t]4r 


< C{d, X)\E\ 2 R 2 £ 2 r . 

The combination of both r-ranges yields 


(27) - 


\ 1/2 


|V£g - V(b%)i(xi + &)\ 2 dx < C(d, A)|£| min{l, 2 R/r}e 2R . 


In total, we see that 


^E V/E + £e — (b E )i{ x i + <Pi) 


is the desired function (note in particular that the last term is a-harmonic), provided 
we choose C\ to be the constant appearing in (1271) . □ 


We now establish existence of second-order correctors by means of the previous 
lemma. 


Proof of Theorem [H We just need to construct an “initial” second-order corrector 
ip r E subject to the properties of Lemma [HI then Lemma |T3] yields a sequence 
(ip E r °) m which is a Cauchy sequence in H 1 {B R ) for every R > 0 due to the last 
estimate in the lemma and our assumption which implies summability of £ 2 m r 0 - 
Thus, the limit tp E satisfies the equation m in the whole space, depends linearly 
on E, and satisfies the estimate 

\ 1/2 oo 

\Vip E \ 2 dx) < Ci(d,X)\E\ ^2 min{l,2 7n ro/r}e 2mr . 0 

' m =0 

for any r > ro. 

To construct ip 7 ^, just use Lax-Milgram to find the solution ip 7 ^ on R d with 
square-integrable gradient (unique up to an additive constant) to the equation 

- V- aVip^ = Eij V-[xs ro (o’ij + Oji + a(<piCj + (pjef))}. 
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Obviously, after fixing the additive constant appropriately tp r B depends linearly on 
E. Furthermore, we have the energy estimate 

[ |VV £| 2 dx < C(X)\E\ 2 [ | XSro <7| 2 + IxB^acj)] 2 dx, 

J R d J R d 

i.e. for any r > ro 

f \Vr E °\ 2 dx<C(d,X)\E\ 2 [ \^\ 2 + \a\ 2 dx 

J B r J B rQ 

and therefore 

/ IVV4 °| 2 dx < C(d,X)\E\ 2 r- d s 2 ro r 2 + d 

J B r 

< C{d , X)\E\ 2 r 2 min{l, (r 0 /r) 2 }e 2 o . 

We note that this provides the starting point for Lemma fldl possibly after enlarging 
the constant C\ in the statement thereof. □ 


3.4. Proof of the C 2,a Liouville principle. The C 2,a Liouville principle (Corol- 
lary[Hl) is an easy consequence of our large-scale excess-decay estimate (Theorem[7]). 


Proof of Corollary{8i Let a G (0,1) be such that 


1 


lim „„ , 

R—¥oo 


lul" dx 


Br 


1/2 


= 0 


holds. By the Caccioppoli estimate, we deduce 


1 


lim 

R —>oo 


/ IV' 

J Br 


u\ 2 dx 


1/2 


= 0 . 


Fix r > ro- The excess-decay estimate from Theorem [3 yields together with the 
trivial bound Exc 2 (f?) < f BR |Vu | 2 dx that 

/ r \ 2+2a 

Exc 2 (r) < C(d, A, a) J Exc 2 (i?) 


< 


C(d, X, a)r 2+2a \Vu\ 2 dx^j ' j . 


Passing to the limit R —> oo, we deduce that 

Exc 2 (r) = 0 


holds for every r > ro- Therefore, on every B r with r > ro, V - u can be represented 
exactly as the derivative of a corrected polynomial of second order (since the infi- 
mum in the definition of Exc 2 is actually attained, as noted at the beginning of the 
proof of Lemma flT)l) . i.e. we have 

V« = V&f (xi + <f>i) + VElj(xiXj + Xi4>j + < j>iXj + ipij) 

in B r for some b r G and some E r G £. It is not difficult to show that for r 
large enough, the b r and E r are actually independent of r and define some common 
b e l d and E G £: For example, one may use Lemma [9] to compare the 6 r , E r 
for two different radii rr,r 2 > ro; the estimate for \b ri — b 1 ' 2 \ and | E ri — E r2 \ then 
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contains the factor Exc 2 (max(ri, r 2 )) and is therefore zero. Moreover, the gradient 
Vu determines the function u itself up to a constant, i.e. we have 

u = a + bi{xi + fa) + Eij(xiXj + Xi(j)j + fa Xj + fa]) 

for some ogK, some b £ K d , and some E £ £ C M. dxd . □ 


4. A C k ’ a Large-Scale Regularity Theory for Elliptic Equations 
with Random Coefficients 

We now generalize our proofs from the C 2, “ case in order to correct polynomi¬ 
als of order k and obtain our C k,a large-scale regularity theory. We proceed by 
induction in k. 

In order to establish our C k,a regularity theory, let us first show Proposition [2] 
which - like the proof of Proposition [ 6 ] in the C 2,a case - only requires a simple 
computation. 


Proof of Proposition [H Making use of the fact that we have ( ahom)ijdidjP = 0 (in 
the third step below), we obtain 

- V • {a t VdiP) 

= (V • o-i) • VfaP 

= Qi ■ Vd,; P 

= afai + Vfa) ■ VfaP 

iv-(a,Pa( e , + v^). 


This yields 


V ■ ((faa - (Ji)S7diP) 

= V • a(faVdiP + d l Pe l + diPVfa ) 

= V ■ aV(P + fadiP), 

which together with (j9j) implies our proposition. □ 

4.1. The C k,a excess-decay estimate. To establish our C k,a excess-decay esti¬ 
mate, we make use of the following lemma, which essentially generalizes Theorem [3] 
to correctors that are only available on balls Bn. 

Lemma 14. Let d > 2 and k >2. Suppose that Theorem]]] holds for orders 2, ..., 
k — 1, and set tpp = 0 for first-order polynomials P to simplify notation. For any 
P £ 'Pa hom > denote by ifp a solution to the equation © on the ball Bp (without 
boundary conditions); assume that the ifjp depend linearly on P. Set 

( 28 ) £ J,rR~ SU P P~ (k ~ 1) ( max -f |V^| 2 dx 
r<p<R \ p ZV* hom ,\\P\\=lJ Bp 
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For an a-harmonic function u in Br, consider the kth-order excess 
(29) 


Exc fc (r) := 


inf 

P K £P, 


~ k— 1 

f j- Vn — V( ^~^(-P k + <t>idiP K + V’Pk) + (-Pfc + fiidiPk + TpPtS) 

a hom'*Pr ^ — '| 


dec. 


For any 0 < a < 1 there exists a constant £ m in > 0 depending only on d, k, A, 
and a such that the following assertion holds: 

Suppose that ro > 0 satisfies £2 , ro + £ f, ro r < £ min ■ Then for all r £ [ro, R] the 
C k,a excess-decay estimate 

-- / 7* \ 2 {k-l)-\-2ot -- 

(30) Excfc(r) < C(d, k, A, a) y—J Exc k(R) 

is satisfied. 

Note that the infimum in (1291) is actually attained, as the average integral in 
the definition of Exc 2 (p) is a quadratic functional of P K . Denote by pp> mm a 
corresponding optimal choice of P K in (USD. We then have the estimates 


k 

(31) R 2 ^~ l) \\P r F min - p?’ min II 2 < C[d, k, A, a)Sre k {R) 

K=1 


and 
(32) 

Proof of Theorem [3 Once we have shown Theorem [TJ Theorem [3] obviously follows 
from Lemma 1141 bv setting ifp k := ifp k , with ipp k being the fcth-order corrector 
whose existence is established in Theorem [T| □ 


k „ 

V R 2 ^- 1 ') | |P£’ min 11 2 < C(d, k, A, a)j- |Vn| 2 dx. 
K =i Jb * 


The following lemma is essentially a special case of our C fc,Q large-scale excess- 
decay estimate Lemma [14J it entails the general case of Lemma Qj] cf. below. 


Lemma 15. Let d > 2, k > 2, and let R, r > 0 satisfy r < R/4, and £ 2 ,r < 
£o(d, k — 1, A), with £o(d, k — 1, A) being the constant from Theorem\7\for the orders 
2, ..., k — 1. Assume that Theorem] 7] holds for orders 2, ..., k — 1, and let ipp = 0 
for linear polynomials P in order to simplify notation. For any P £ Pa hom > denote 
by i[>p a solution to the equation © on the ball Br (without boundary conditions); 
assume that ipp depends linearly on P. For an a-harmonic function u on Br, 
consider again again the kth-order excess (ESI). Then the excess on the smaller ball 
B r is estimated in terms of the excess on the larger ball Br and our quantities £ 2 ,r 
and Vfjp: We have 


Excfc(r) < C(d,k, X)F,xck(R) 


(5) 


2k 


J 2 /(d+ 1) 2 

" 2 ,R 


+ ir 2(fe_1) max 

PePl, ,l|P|l=l.j 


f B N*p \ 2 A Q' 


Before proving Lemma fTKl we would like to show how it implies Lemma fill 
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Proof of Lemma [7/] First choose 0 < 9 < 1/4 so small that the strict inequality 
C(d , k, A )9 2k < g 2 t k ~ 1 )+ 2a i s satisfied (with C(d, k, A) being the constant from 
Lemma fl5l) . Then, choose the threshold s m i n for C 2 ,r 0 + ro R so sma, U that the 
estimate 


C{d,k, A) 


9 2k + (4[ro +1? + e2 


il),ro,R 


^ Q 2 {k— l)-\- 2 a 


holds. 

Let M be the largest integer for which d M R > r holds. Applying Lemma [T5l 
inductively with R m := 6 m ~ l R , r m := 9 m R for 1 < m < M, we infer 

Exc k (9 M R) < ( 0 2 (fc - 1 )+ 2 “) M Ekc fe (l?). 


Since we have trivially 


Excfc(r) < 


- 

\tm 


—d 


Exc k(r M ) 


and since by definition of M we have r > 9 vm and thus 9 AI <9 1 (where we 
recall 9 = 6 (d, k, A, a)), we infer 


Exc fc (r) < C(d, k, A, a) 


\ 2(fc—l)+2o; 


ExCfc(P). 


It remains to show the estimates for || pr,mm _ pR,mm || ag wed ag ^hg F> ounc l s for 

|pr,mm|| Tq do g0j ^ ug firgt estimate t he differences \\pR^,min _ pr m ,rnin || of 

two successive polynomials. We have the estimate 


k -1 


/ v { p * m ,mm - p * m ,min +&’ 
• ' B r ~~ - -1 

4>idi(P k 


ram _ ppm ,rmn^ _|_ ^ 


° r m K — 1 


0 


VP, 


Rm,min r> r m ,min 
Ji„ 


k k 

k -1 


- P r 


i\ , 7 

) +Vh 


*-pr 


’_p r 

r L. 


< 2 / Vn - V ^ (p;-’ mi " + fadiP r K ™' min + 


K =1 


-VP, 


0 


-tV-m .min 


1 r\ 73)r m ,razn 

• ('piOir f, -|“ ’IfJ pr-m ,min 


'V'f 


da; 


/. fe— 1 

+ 2 / Vll - V V (pK^rnin + ^.Q.pRn^min + 


k=1 


-VP 


0 




d 


< 2Exc fc (r m ) + 2 

C(d,fc, A,a)(^) 

I'm \ 2 ( fc_1 ) 


< 


R 
r r 

Tm \ 2 ( fe — 1)+ 2q: 


Exc/- {R rn ) 


+ Ip pRm ,min 


Rn 


Exc fe (P) + C(d, fc, A, a) 6 »" d 

R 


dx 


2 (fc-l)+2a 


ExCfe(P) 


< 


( r \ Z[K— 1 ) 

^f) (d 2 “) m Exc fc (P). 
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From Lemma [THl below, we thus obtain 


^RK-l\\pRrn,min ^ pr m ,min || < (7((i, fc, A, tt) (#“ ) m y ExCfc (P). 

K = 1 

A similar estimate for the last increment Yl,=i — pr,mm || can j-,e 

derived analogously. Taking the sum with respect to m and recalling that R\ = R 
and r m = R m+ i, we finally deduce 


^ ^ pn—1 11 pR,min _ pr,min 11 

K.— 1 

M r= - 

< C(d, k, A, a) Y ( 0a ) m \/Exc fc (i?) 


< C(d, k, A, a)y ExCfc(P). 


It only remains to establish the last estimate for ||P£’ mm ||. By the previous esti¬ 
mate, it is sufficient to prove the corresponding bound for ||P^ ,mm ||. This in turn 
is a consequence of the obvious inequality 


fc-i 


/ v E ( p ?’ min +<t>i d i p ?' min + v>p“) 

JBr K=1 V 

, V7 ( T-)R,min . / o T-iR.min . 7 A 

+ V (Pfc + <l>idiP k + ‘IppR.min J 

< 2ExCfe(P) + 2P |Vu | 2 dx < 4 J- |Vu | 2 da; 
JBr J Br 


dx 


in conjunction with Lemma ITTT1 below. 


□ 


The following lemma quantifies the linear independence of the corrected polyno¬ 
mials P K + 4>idiP K + i\)p K (with 1 < k < fc); it is needed for the previous proof. 

Lemma 16. Suppose that the functions <p and ifp^ (2 < k < k) satisfy 



p -2( K -D 


max 

PeV 2 h om' II P H = 1 



|V0 p| 2 dx < el, 


where £q = £o(d,k) is to be defined in the proof below. Set ipp = 0 for linear 
polynomials P in order to simplify notation. Then for any P K £ Pa hom 0- is K < k) 
we have the estimate 


(33) 


E^ : 

AC=1 


2(k-1)| 


\Pk. 


< 


P k 

C(d, k)j- VV(P, 

Jb p k=1 


fadiPK + ifpJ 


dx. 
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Proof. Poincare’s inequality (with zero mean) and the triangle inequality imply 


' r k 

/J V S (F “ 


+ 4>idiPn + 'ijjp K ) 


1/2 


dx 


> -Jr- - inf 


> 


d/(d) p aCM B 
1 1 


C(d)p 


~ K 

t ^^(Pk, + 4 >idiP K + — a 

K=1 

x 1/2 


2 \V2 

dx 


inf 

aeR 




P K - a 


P K=1 


dx 


/ fc 

■?£, /„ Dw.+fe) 


— a 


2 X 1 / 2 

dec 


On the one hand, by transversality of constant, linear, homogeneous second-order, 
..., and homogeneous fctli-order polynomials we have 


-inf (7 £p k 
P \J Bo 


— a 


1/2 


dx 


> 




On the other hand, we have by the triangle inequality and Poincare’s inequality 
k 

((hi f)i P„ 4- fi) d 1 — n , d.x 


— inf 

p a(ER \JB 


/ ^ {4>id,P K + VtpJ - a 

J B " K=1 
k 


2 \ V 2 

dx 


< C(d, k) 




1/2 


p K MlPJI-r max 

“ P K_1 Peep-,||P||=i y B 


. (£. |v ^ 5 


1/2 


da: 


Putting these estimates together, by boundedness of the integrals in the previous 
line by £q p 2 ( K_1 ) our assertion is established. □ 

Proof of Lemma \15\ In the proof of the lemma, we may assume that 

(34) ExCfc(I?) = / |Vu | 2 dx. 

Jb r 

To see this, recall that the infimum in the definition of ExCfc(P) is actually attained. 
Denote the corresponding choices of P K by P™ n . Replacing u by u — ^2 t Z} 1 (P™ ln + 
(j>id i P™ ln +'ifp™.ird)-{PJf lm +(j)id i PJf lln prnin)), we see that we may indeed assume 
(1M1) : The new function is also a-harmonic due to © and Proposition [2] 

We then apply Lemma [20] below to our function u. This yields an a/j 0m -harmonic 
function Uhom close to u which in particular satisfies 


/ \^u hom \ 2 dx < C(d, X)J- |Vu | 2 dx. 

" B R / 2 J BR 
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By inner regularity theory for elliptic equations with constant coefficients, the cihom- 
harmonic function Uhom satisfies 

k 

(35) |Vtifcom(0)1 +R sup \V 2 u hom \ + Y RK sup |V K+ 1 Uh om | 

Br/4 k= 2 b R /4 

( \ 1/2 

/ |Vu ^ om | 2 c&r ) <C(d,k,X)(f \Vu\ 2 dx 

J b R / 2 J \J Br 

Let pB,Taylor (f or l < n < k) be the term of order k in the Taylor expansion of 
Uhom at To = 0. We now show (for k > 2, as for n = 1 this assertion is trivial) 
that pB,Taylor g V ^ hom . The term-wise Hessian of the Taylor series of Uhom yields 
the Taylor series of V 2 rtfc, om . We now know that a/, om : V 2 tt/ lom = 0; thus, the 
Taylor series of a/, om : V 2 rt/ lom is identically zero and by equating the coefficients 
we deduce ah om ■ y 2 pB^ Ta v lor = o for 2 < n < k. 

As the term-wise derivative of the Taylor series of Uh om yields the Taylor series 
of we obtain by the standard error estimate for the Taylor expansion of 

Vit/tom at Tq = 0 for any x £ B R / 4 the estimate 



< M SU P U ho m\- 


lR/4 


Vu hom (x) - Y VP?’ Tayl ° r (x) 

K,= 1 

Making use of the identity 

k 

(Id+(V</>) t )Vu/ lom - vY( p *’ Tavl ° r + <!>idiP*' Tavlor ) 

K= 1 

k 

+ Y&Vd l pK' Taylor 

k=2 

= (Id+(V</>) 4 ) (vu hom (x) -YVP?’ Taylor (x)) , 


the previous estimate yields in connection with the bound for |V fe+ 1 u^ om | and 
r < i?/4 


r ,. k 

f (Id +{S7cj>) t )Vu hom - V Y(P?’ Taylor + (j>idiP^ Taylor ) 

JB r *=1 

k 2 

+ Y (j) i Vd i PY aylor dx 

< C(d, k, A) ( Y) 2K / |Vu | 2 dx x -f | Id +(V ^)*| 2 dx. 

' PP J Br J B r 

By the Caccioppoli inequality for the a-harmonic function xi+<fii (cf. ©), we have 
(36) -f | Id+(V ^))*| 2 dx < ^ / \x + (fr \ 2 dx < C(d, A)(l + e%). 

JBr r JB.r 
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The approximation property of Uhom + 4>idiUhom in Br /2 from Lemma [20] below 
implies 

j- |Vu - V(uft om + (f>idiU h om)\ 2 dx < C(d, A)e^ /(d+1) £ |Vit | 2 dx. 

Combining the last three estimates and the equality 
k -1 

Vu (p*' Tavlor + 4Hd i pX’ Tavlor + ^ p ^ Tayl \ 


-VP, 


K=1 

^ p R,Taylor + faQ.pR,'Taylor + 
k 

(Id +(V#)V« t « - V £ (p^aylar + frd.pK,Taylor^ 


K—l 


k-1 


+ J2 frVdiP*' Tayl0r - ]T c/ )i Vd i P^ Ta y l0r - J2 V^pR,Taylor - 

k=2 j /t=2 k=2 

+ ^(dhom H” 4 > i^i'^ J hom) H” 0iV diUhom > 


we infer 


» k—l 

/ Vn - V ^ ^pR,Taylor + g.g.pR,Taylor + ^ 


-VP 


0 


/t=l 
R,Taylor 


+ 4>idiP, 


R,Taylor 


i^i± k ' r P£ 

k 


') 


da; 


< 


ijj pR,Taylor 

(Id + (V$*)'Vu hom - V ^ ^pR,Taylor + ^.g.p^Taylor'j 
J B r K=1 

k 2 

+ ^ <j) i Wd i P^ Tavlor dx 

k=2 

C{k)l V |V 2 P^—1 2 |^| 2 dx 
J B r 


k=2 

k—l 




R,Taylor 12 


da; 


< 


+c(k)j- 

S '- re =2 

+ 6f |ViApH ,Taylor | 2 da; 

J B r k 

+ 6 j- |Vit — V(ztft, om H - 4 ) i&i r U'hom) \ dx 

J B r 

+ 6+ |0iV(9i^ om | 2 dx 

JB r 

C(d,k,\)(^Pj (l + £ r)j- |Vm | 2 da; 


+ <7(d, fc) ^ r 2 ^- 1 ) | IP^ TayZor 112 £ 2 
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■C(d,k)^2\\P^ Taylor \\ 2 max ||P||- 2 £ \W dx 
1 max 

p ev* ho 


ac=2 

■ C(d, k)\\P^’ Taylor \\ 2 max ||P|| _2 / |Vi4| 2 dx 

*■_ 

^2 / r \ ~ d r 


- C{d , A)4 /(d+1)5 




+ C{d)r 2 e 2 sup |V 2 M/,, orn | 2 . 

Br/a 

This finally yields in connection with the bounds on S/^Uhom in Br / 4 (cf. (1551) 1 

/ \ 1/2 

which in particular imply ||P^ To ^ or || < C(d,k, A)P 1 " K U Br |Vw | 2 da;J 
/» fe —1 

f Vu-VV (pR,Taylor + ^.g.pR,Taylor + ^ pRTayl \ 

J Br K= 1 V 7 

- V (P fe H ’ Ta ^° r + <\> i d r P*' TaVl ° r + dx 

<C(d,k, A) (l + £ l)J- l Vu | 2 dx 

+ C(d,fc,A)e 2 ^(L) ” 1 \Wu\ 2 dx 

k=2 ' ■'•Bh 

« fc— 1 « 

+ C(d, fc, A)/ |Vm | 2 dx y P 2 ^ : ~ 1 ' max ||P||~ 2 X IVi/ApI 2 da; 

■'Pr k =2 PeP K 

+ C(d,k,X)-f |Vm | 2 dx x p- 2 ^ -1 ) max ||P||~ 2 / |Vi/ip | 2 dx 
Jb r p ^ h om J Br 

+ c(d,A) 4 /(d+ 1 ) 2 (L)^£ | Vl 

+ C(d,X)r 2 £ 2 r R~ 2 -f |Vit | 2 dx 

Jb r 

r [ , r \ 2k 

< C(d, k, A)-/- |Vw | 2 dx (- 

JB» L v -^ 7 


7u| 2 da; 


+ £ 2 7 p +1)2 + P- 2(fc - 1} max 

V P 6 ^\ om -H P il= 1 -/PR 




where in the last step we have used the inequality £ 2 < (-y) d £p < ( p ) d £p^ +1) 
and £p < £ 2 : p as well as (flOl) for 2 < k < k — 1. Our new estimate now implies the 
desired bound. □ 


4.2. The C k 1,1 excess-decay estimate. Like in the C 2 ' a case, we now show 
how the C k ’ a excess-decay estimate for the fcth-order excess Excfc (cf. Lemma fldl) 
entails a C k ~ 1,x excess-decay estimate for the (k — l)th-order excess Excfc_i. 

Lemma 17. Let d > 2, k > 2, and R > 0. Assume that Theorem[ 7] holds for the 
orders 2, ..., k — 1, and let ifp = 0 for linear polynomials P in order to simplify 






26 


JULIAN FISCHER AND FELIX OTTO 


notation. For any P £ 'Pa hom > denote by ipp a solution to the equation © on the 
ball Br (without boundary conditions); assume that the ifp depend linearly on P. 
Then there exists a constant £ m in > 0 depending only on d, k, and A such that the 
following assertion holds: 

Suppose ro £ (0, R] is so large that £ 2 ,r 0 < £ m in and 


sup p 
r 0 <p<R 


-(fc-i) 


max 

Pert ,||P||=i 


L 


1/2 


|V-0p| 2 dx 


_ '-■min 


hold. Let u be an a-harmonic function on Bp. Then there exist Pff £ P” h 
(1 < k < k — 1) for which the estimate 


L 


k -1 

S/U-VJ2 ( p f + MiP? + 4>p*) 

K=1 

/ r\ 2 (M f 

C(d,k,X ){-) f |Vu| 


dx 


dx 


Br 


holds for any r £ [tq , R] ■ Furthermore, the P^ depend linearly on u and satisfy 


k—1 

^i? 2 ( K -l)||pR ||2 < C (d,k,\)l |Vm| 2 dx. 

J Br 


Proof. In Lemma [Til fix a := 1/2. We then easily verify that Lemma [T-H is applica¬ 
ble in our situation. Set Pff := pr 0 ,m»n. ^his implies that the Pf 1 depend linearly 
on u. The estimate (15^1) takes the form 


k 

E 

K.— 1 


R 2 ( K -i)^ p R ^2 < C (d,k,X)f |Vm| 2 dx. 

JBr 


Furthermore, applying Lemma [IT] with ro playing the role of r and r playing the 
role of R, we deduce from m 


5>2(«—D||P«_P^|| 2 


< C(d, k, A)ExCfe(r) 

r X 2(fe-l)+2 

R 

C(d,k,X )(^ 


m / v \ 

< C{d,k, A)(-J 


< 


\ 2(/c—1) —(— 2 ot 

R.) 


Exc k (R) 



|Vm| 2 dx. 
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We now estimate 

fc-i 


/ Vu _ V £ (P* + fc^P* + ^) 

■ /B r * =1 

« fc—1 

< 2>j- Vti - V 53 (P” + fadiPl'™™ + l/jpr.min ) 

2 

V / nr,min . / nr,min , 7 \ 

(Pfc + <^ P fc + l/jpT.minj 

+ 3 f v(p” + ^P” + 2 

+ 3 £ v E ( p : ,min - p *+^( p « ,min - p *)+^ 


dx 


dx 




dx 


< 3ExCfe(r) 


+ C(d,fc)||P fe r ’ mirl || 2 r 2 ( fc - 2 )^ |^| 2 +r 2 |Id+(V</>)T dx 

+ r -2(k-2) max f |y^ p | 2 da-] 

pe 7 ’ih» m >ii j> ii= i ^ P y 

fc -1 „ 

+ C(d, fc) 53 r 2< - K_1 ^||P^’ mi " — Pk\\ 2 t | Id + (V0)*| 2 dx 

JB t 

|| 2 max / r 2 ( K “ 2 )u| 2 _|_ |VV’p| 2 dx 


K=1 

fc-1 


+ C(d, k ) V | |PP™" - P^H 2 max 

v ' ^ " K K " Per? ,||P||=iJ B 

tt=2 a hom ’ 11 

J10I30I361 

C(d, k, A) J Excfc(P) 

+ C(d, fc, A) | |pP min 11 V (fc_1) (e 2 + (1 + s 2 r ) + ej ^ 

+ C(d , fc, A) ^ r 2(«-i) | |pr,min _ pP||2 (1 + ^ 

K = 1 

+ C-(d, fc, A) 5] r 2(«-D ||P~ _ pP||2 (e 2 + e 2 r) 


< 


k=2 

/ r \ 2(fc—l)+2a __ . 

C(d,fc,A) (-) Excfc(P) + C(d, k, A)||P^’ mm || 2 r 2 ^ fe_1 ^ 


,P 

fc-1 


+ C(d, fc, A) 53 r2(K_1) | |P^ min - P^ 112 


/■C=l 


In conjunction with the two previous estimates, we infer 

fc -1 2 

dx 

K,= 1 


„ fc-i 

f vu - v 53 (p«++v> P «) 

J B r k-—i 


< C(d, k, A) 


/ T* \ 2(fc — l)+2o: X / T \ 2(fc—1) / y \ 2(fc—1)-|-2 q'\ / T* \ —1) - | - 2q! 

(p) + ((p) + (p) ) + (p) 
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X 



dx. 


Our lemma is therefore established. 


□ 


4.3. Construction of correctors of order fc. Using the C k ~ 1,1 theory estab¬ 
lished in the previous subsection, we now proceed to the construction of our fcth- 
order corrector. The following lemma provides the inductive step; starting from a 
function which acts as a fcth-order corrector on a ball Br, we construct a function 
acting as a fcth-order corrector on the ball B 2 r. 


Lemma 18. Let d > 2, fc > 2, and assume that Theorem\7\ holds for the orders 2, 
.... fc — 1. Let r 0 > 0 satisfy the estimate e 2 , ro fs £o> where £o = £o(d, fc, A) is to be 
chosen in the proof below. Then the following implication holds: 

Let R = 2 M ro for some M £ No- Suppose that for every P £ V k we have a 
solution if>p to the equation 

— V- aVipp = V ixBni&a - cq)V diP) 
subject to the growth condition 

s 1/2 M 

\Vipp\ 2 dx J < Ci(d, fc, A)||P|| ^2 min{l,2 rn ro/r}e 2 m T . 0 

' m —0 

for all r > ro, where Ci(d, fc, A) is a sufficiently large constant to be chosen in the 
proof below. Assume furthermore that ip R depends linearly on P. 

Then for every P £ V k there exists a solution ipjf 1 to the equation 

- V- aVifp R = V-(xB 2R (<t>ia - aJ'VdiP) 



subject to the growth condition 

\ 1/2 


r -{k- 1 ) 


/ r \ 1/2 M -\-1 

(f \V^ 2 P R \ 2 dx) <C 1 (d,fc,A)||P|| V min{l,2 m r 0 /r}e 2 

VB r / m=0 


for all r >vq. Furthermore, , f) 2R depends linearly on P and we have 
r_(fc_1) (/ “ V ^! 2 dx ) 7 < C'i(d,fc,A)||P||e 2 M+ 


Proof. To establish the lemma, we first note that the assumptions of the lemma 
ensure that the C* -1,1 excess-decay lemma (Lemma flTl) is applicable on Br with 
i[)p := if R . To see this, we estimate for any r £ [ro,i?] 


r -(fc-1) 



< C 1 {d,k,X)\\P\\£ 2 ,r 0 


< Ci(d, fc, A)||P||e 0 . 


By choosing £ 0 > 0 small enough depending only on d, fc, A, and C\ (which is to be 
chosen at the end of this proof), we can ensure that the assumption of Lemma [T71 
regarding smallness of e^ ro R is satisfied. 

We now turn to the construction of if™ — ipp and to that purpose denote by 
the weak solution on with zero mean in B 2 r and square integrable gradi¬ 
ent, whose existence and uniqueness follows by the Lax-Milgram theorem, to the 
problem 


V- aV£p = V-(xB 2R -B R {(pia - a^S/diP). 
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Obviously, £p depends linearly on P. Furthermore, by ellipticity we have the 
estimate 

ive£i 2 dx 

< C(d, A) sup |V 2 P| ( [ Xb 2 R -b r (.\M 2 + W\ 2 ) dx) ( [ |V^| 2 dx) 

B 2R \JR d / \JR d / 


which gives 


( [ |V£p| 2 dx) < C(d, A) sup |V 2 P| f ( \cp\ 2 + |tr| 2 dx 

\jR d ) B 2R \J B 2R 


\ 1/2 


The last estimate in turn implies 
(37) 


[ |V^| 2 dx < C{d,k,\)\\P\\ 2 R 2{k - 2) £ 2 2R R 2+d . 

J R d 

We now obtain ip 2R — ipp by modifying £p by an a-harmonic function of degree 
k — 1. As £p is a-harmonic in Bp, Lemma 1171 now implies the existence of some 
P R P £ V K for 1 < k < k — 1 which depend linearly on P and for which the estimates 


(38) ||P* F || 2 < C{d,k,X)R~ 2 ^l |V£p| 2 dx V C(d,k,X)\\P\\ 2 R 2 ^e 2 R 

Jb r 


and 


» k— 1 

f - V £ (p R p + oApKp + ,p p n p ) 

J B r _i 


K=1 

/ r \ 2 M r 

< C(d,k, A)(-j ^ |VC 


dx 


R ' 2 dx 


J33 


< C(d,k,X)\\P\\ 2 r 2 ^e 2 2R 


dx 


hold for all r G [ro,P]. 

Furthermore, we have for r > R 

« k—1 

i (p r p +<P idiP R P + ip PRf 

■' B r re=l 

(m / r 

< C(d,k,\)lr~ d |Vf^| 2 di + ||P'y|| ! (l + 4) 

+ y),»-')||p« P ii 2 (i+ E 2 r + 4 r )) 

J 

C(d,k,X)\\P\\ 2 R 2 ^ + l+e 2 r +(l+S2r+S 2 , r )(^y (k ^ 

< C{d,k,X)\\P\\ 2 r 2 ^R 2 e 2 R . 


e 2 R 
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The combination of both r-ranges yields 



< C(d, k, X)\\P\\mm{l,2R/r}e2R. 
In total, we see that 


k —1 

V’p* := i’p + £p — + V’pSp) 

K=1 

is the desired function (note in particular that the last term is o-harmonic), provided 
we choose C\ to be the constant appearing in (15(71) . □ 


We now establish existence of fcth-order correctors by means of the previous 
lemma. 


Proof of TheoremUl We just need to construct an “initial” /cth-order corrector ip'p 
subject to the properties of Lemma 1T51 then Lemma [lH] yields a sequence (ifp r °) m 
which (after subtracting appropriate constants) is a Cauchy sequence in PI 1 {Bp) 
for every R > 0 due to the last estimate in the lemma and our assumption (Hj) 
which implies summability of e 2 m ro . Thus, the limit ifp satisfies the equation © 
in the whole space, depends linearly on P, and satisfies the estimate 

\ 1/2 oo 

\Vipp\ 2 dx) < Ci{d,k, A)||P|| ^ min{l,2 m r 0 /r}e 2 ">r- 0 

' m =0 

<Ci(d,M)||P||e 2 , r 



for any r > tq. 

To construct i/’p , we use Lax-Milgram to find the (unique) solution ip r p on M d 
with square-integrable gradient and zero mean on B ro to the equation 


- V- aVi/>p° = V-( XBro ((fid - <Ti)VdiP). 

Obviously, iffS depends linearly on P. Furthermore, we have the energy estimate 
[ IVi/jp) 2 dx 

< C(d, A) sup |V 2 P| ( [ \XB ro a(t)\ 2 + \xB ro <j\ 2 dx) ([ |V0p°| 2 dx) 

B ro \jRd J \jRd J 

We therefore get 

( [ IV^p 0 ) 2 dx^\ < C(d, A) sup |V 2 P| ( f \4>\ 2 + \a\ 2 dx J 
\JR d J B r(J \JB ro J 

This yields in particular for any r > ro 

/ \V^ r P °\ 2 dx<C(d,k,X)\\P\\ 2 r 2{k ~ 2) [ \0\ 2 + \<j\ 2 dx 

J J 


1/2 
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and therefore 

/ W| 2 dx < C(d,k,\)\\P\\ 2 r- d rf k ~ 2) el o rl +d 
J B r 

<C(d,k,X)\\P\\ 2 r 2 ^mm{l,(r 0 /r) 2 }e 2 ro . 

We note that this provides the starting point for Lemma[T8l possibly after enlarging 
the constant C\ in the statement thereof. □ 


4.4. Proof of the fcth-order Liouville principle. Like in the C 2,a case, the 
C k ’ a Liouville principle (Lemma ITTFlbelowl is an easy consequence of our large-scale 
excess-decay estimate (Theorem[5]). The fcth-order Liouville principle (Corollary 0]) 
in turn is an easy consequence of the C k+1 ’ a Liouville principle. 

Lemma 19. Let d > 2, k > 2, and suppose that the assumption © is satisfied. 
Then the following property holds: Any a-harmonic function u satisfying the growth 
condition 

( 4 °) fiE+L (£ M* dx ) = 0 

for some a £ (0,1) is of the form 

k 

U = a + bi(xi + 4>i) + y^(P K + (t>idiP K + ipp K ) 

k—2 

with some a € R, b £ R d , and P K £ 'Pa hom f or 2 < n < k (i.e. P K is a homogeneous 
ahom-harmonic polynomial of degree k). Here, the %f)p denote the higher-order cor¬ 
rectors whose existence is guaranteed by Theorem[ T[ 

Proof of Corollary [7} Obviously, (fl3l) entails (HU1) with k + 1 in place of k and e.g. 
a\=\. By Lemma fT9l any a-harmonic function u subject to condition (fTSl) must 
be of the form 

k +1 

( 41 ) u = a + bfixi + fii) + ^(P K + fiidiP^ + 4>p k ) 

k—2 

with some sgK, b £ R d , and P K £ Va hom for 2 < k < k + 1. Our stronger growth 
condition m however shows that we have Pkp \ = Since the ci, grow sublinearly, 
cf. ©, and since ifp k+1 grows slower than a polynomial of degree k + 1, cf. (fTUl) . 
we see that for large |x| the term Pk+i would be the dominating term in (14T1) if it 
were nonzero, contradicting our growth condition (1131) . □ 


Proof of LemmaUM Let a £ (0,1) be such that 


1 


liminf . , 
R—*■ oo R k +° 


'/ 

J Br 


|m| dx 


1/2 


= 0 


holds. By the Caccioppoli estimate, we deduce 


1 


liminf . , , 

R-s-oo 


j- |Vu| 2 dx ^ 


1/2 


= 0. 
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Fix r > tq. The excess-decay estimate from Theorem [3] together with the trivial 
bound Exc fc (P) < f g |Vu| 2 dx yields 


ExCfc(r) < C(d,k,X 1 a) 


/ 'T \ 2(/c—1)+2 ck 

vpJ 


Exc fe (i?) 


< C(d,k,X,a)r 2{k - 1)+2 ° 


R k ~ l+a (1 |Vm| 


Passing to the lim inf R —> oo, we deduce that 

Excfc(r) = 0 



holds for every r > r 0 . Therefore, on every B r with r > rg, Vu can be represented 
exactly as the derivative of a corrected polynomial of A'tli order (since the infimum 
in the definition of ExCfc is actually attained, as noted at the beginning of the proof 
of Lemma 1 151) . i.e. we have 


k 

Vu = Vb\{ Xi + fa) + Vj2( pr . + 

k=2 

in B r for some b r £ M. d and some P£ £ P a hom (2 < n < k)\ recall that we have 
used the convention ipp = 0 for linear polynomials P. It is not difficult to show 
that for r large enough, the b r and P£ are actually independent of r and define 
some common b £ I d and P K £ P a hom - For example, one may use Lemma 1X71 to 
compare the b r , P£ for two different radii rr,r 2 > ro; the estimate for \b ri — b r2 \ 
and ||P^ 1 — P^ 2 1| then contains the factor Excjfc(max(ri, r 2 )) and is therefore zero. 
Moreover, the gradient Vu determines the function u itself up to a constant, i.e. 
we have 

k 

U = a + bi (xi + <j>i) + '^2 {Pk+ (t>idiP K + i/jpJ 

K ,=2 


for some a £ R, b £ l d , and P K £ p a hom (2 < k < k). 


□ 


Appendix A. Approximation of ci-Harmonic Functions by Corrected 

niiom - H armonic Functions 

Our proofs make use of the following lemma, which is implicitly derived in the 
course of the proof of Lemma 2 in (7j. For the reader’s convenience, we recall its 
proof here. 

The lemma essentially states that an a-harmonic function u on a ball Bp may 
be approximated on the ball Bp / 2 up to a small error (of order £^ rf+1 ^ ) by an 
appropriate a/jom-harmonic function Uhom and correcting this function Uhom using 
the first-order corrector fa. 

The purpose of the lemma is the same as in classical elliptic regularity theory: 
The function Uhom satisfies an elliptic equation with constant coefficients, i.e. it 
is smooth and good estimates for its higher derivatives are available. In our proof 
above, we show by means of the present lemma that this high regularity of Uhom 
transfers (in an appropriate sense) to u itself. 
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Lemma 20 . Let R > 0 and let u be a-harmonic on Br. Suppose that £r < 1 (with 
er as defined in d3j ). Then there exists an ahom-harmonic function iihom on Br/ 2 
satisfying the following two properties: First, we have the energy estimate 


(42) 


f I'Vuhoml' 2 dx < C{d, X)J- |Vu | 2 dx. 
JB r/2 J b r 


Second, the “corrected” function Uho 
the sense that 


■ (fidiUhom is a good approximation for u in 


L 


+ 4>idiUhom)\ dx ^ C(d,X)£ji 


2 /(d+l) 


R/2 


■/ 

Jb r 


|Vu | 2 dx. 


Proof. Choose some R' G [jR, R] for which 


(43) 


R'-f \Vu\ 2 dS < C(d)l |Vu | 2 dx 

JdB R , Jb r 


holds. Let Uhom be the a/jom-harmonic function in Br , which coincides with u on 
dB r'. Testing the equation — V- ahom^Uhom = 0 with Uhom. — u (note that this test 
function is admissible since we have Uhom — u = 0 on 8Br> ), we infer by ellipticity 
of a and (in the second step) Young’s inequality 


/ |Vu ftom | 2 dx < C(X)f \Vu\\Vu hom \ dx 
i Bn/ J B Tf/ 


(44) 


< 




\Vu horn \ 2 dx + C(A)/ |Vii | 2 dx 


J E 


R' 


which because of R/2 < R’ < R gives the desired energy estimate. It remains to 
establish the approximation property of Uhom + 4>idiUhom- 

Denote by go : R —>■ K a smooth function with r/o(s) = 1 for s > 1 and r/o(s) = 0 
for s < 0. Let 0 < p < R/ 4 and set p(x) := po(2(R' — p/2 — \x\)/p). Note that we 
have |V? 7 | < C{d)/p as well as p = 0 outside of B R i_ p / 2 and rj = 1 in Br>- p . Due 
to p < R/4, we also have R' — p > R/2. We will optimize in this “boundary layer 
thickness” p at the end of the proof. 

Let us abbreviate 


V ’.— U Uhom r](fidiUhom • 

where the purpose of ?y is to have v = 0 on 8Br’. The desired approximation 
property of Uhom + (fidiUhom as stated in the lemma will be a consequence of an 
appropriate energy estimate for v (recall that we have 77 = 1 in Br/ 2 since p < R/4 
and R’ > 3.R/4). 

To derive this energy estimate, we would like to show that v is approximately 
a-harmonic. We first compute using the fact that u and Xi + (/>i are a-harmonic (cf. 
©) 

— V- aVv 

= - V- oVu + V-(l - rf)aS7uhom + V- a(ei + V fi^pdiUhom + V- 0,aV (pdiU h om) 

= V-(l - r])aX7u hom + a(ej + V<fo) • V (qdiU horn ) + V- (fiaS7 (qdiU hom ) 

= ^*(1 Ohom)^Uhom T (^(c^ -f- V (f>j ) R/iom^z) ‘ ^( pdiUhom ) 

+ V • c^aV (pdiU h om) 5 
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where in the last step we have used the ah, om -harmonicity of Uhom in form of the 
equality - V-(l - r])ahom'Vuhom - ah om ei ■ V {rjdiUhom) = 0. Taking into account 
the formula a(ej + V0i) — ahomei = V- cq (cf. Q),©) and the fact that 

(V • <ji) • Vw = dkPijkdjW = dk(&ijkdjw) = - d k (a ik jdjW ) = - V-(a z \/w) 

holds for any function w by skew-symmetry of cq, we may rewrite the right-hand 
side in divergence form: 

- V- aVv = V-(l - rf)(a - ahom)Vuhom + V-(<&a - <Ti)V(rjdiUhom)■ 

Testing the weak formulation of this equation with v (recall that v = 0 on dB ri) 
and using the ellipticity of a, we deduce using Young’s inequality and the properties 
of p 


|Vu| 2 dx 


b K 


^ C(A) / |(1 ^?)(a Q*hom)^dlhom\ d" 1<7? | \ ^7 {jjdiUhom) |~ dx 

Jb r , 

< C(d, A) f |1 - r]\ 2 \Vu horn \ 2 dx 

Jb r , 

+ C{d, A) f {\(p\ 2 + \a\ 2 )(\S/r]\ 2 \^/u horn \ 2 + T] 2 \\7 2 u h om\ 2 ) dx 
Jb r , 

< C{d, A) [ | V Uhom | ^ dx 

d Br/ —B R /_ p 

+ C(d, A) sup iVu^oml 2 + IV^oml 2 ) [ \(/)\ 2 + \a\ 2 dx. 

Br'-p /2 \P / Jb r i 


Since our function Uhom is ah om -harmonic, we have the regularity estimates 

sup ( \\S7u hom \ 2 + \V 2 u hom \ 2 ) < C ( d ^ sup / \S7u hom \ 2 dx, 
B R ’- p /2 \P / P V^B r i_ p / 2 JB p / 2 {y) 

([ |Vu fcom | p dz) <C(d, X) f \V tan u hom \ 2 dS, 

\Jb r , J JdB R , 

where p := 2 d/{d — 1): The first estimate is a standard constant-coefficient in¬ 
terior regularity estimate (which is a consequence e.g. of an iterative applica¬ 
tion of Theorem 4.9 in [6] and the Sobolev embedding). The second estimate 
follows by combining 1) the existence of an extension u of Uhom subject to the 
estimate ||Vu|| iP ( B ,) < C{d)\\S7 tan Uhom\\i, 2 (dB R ,) and 2) the Calderon-Zygmund 
estimate on Bri , which reads ||Viy||iP(B ,) < C(d, A)||Vu||lp(b r /) for any solution 
w £ H 1 (Br') with w — u £ Hq(Bri) to the equation — V - ahom^w = 0. For the 
latter estimate, see Theorem 7.1 in [8]. 

Using these regularity estimates, the equality V tan tt/ lom = V tan u on 8Br, as 
well as the obvious inequality 


sup 

v^b r i _ P i 2 


£ 


b p/2(v) 


dx £ 



dx , 
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we infer by p < R' /4 and 3i?/4 < R' < R 



I Vu| 2 dx < C(d, A)| B r , - B w _ p I 1 ” 2 /** ( f \Vu hom \ p dx ) 

\Jb r ,-B r ,_ p ) 

dx ■ {R') d l 101 
J Bp, 


1 / f?' \ f 

+ C(d,A)-^-) f |Vn fto 


Jb h , 


2 + |cr| 2 dx 


? C(d,X)p 1/d R' {d 1)/d [ \V tan u\ 2 dS 
JdB R , 

( jdI \ d+2 

—J J |Vm| 2 dx 

GU / n \l/d r 

< C(d, A) (—) J b |Vw| 2 dx 

( TDf \ d+2 „ 

—J J |Vit| 2 dx. 

We optimize in p by choosing p := n' (which thanks to the assumption 

£r < 1 is admissible in the sense of p < \R'). This yields 



|Vw| 2 dx < C{d, A)4 /(d+1)2 



which together with the estimate 3i?/4 < R' < R and rj = 1 in B R / 2 proves the 
desired approximation result. □ 


Appendix B. Failure of Liouville principle for smooth uniformly 

ELLIPTIC COEFFICIENT FIELDS 

We now provide the argument that smoothness of a uniformly elliptic coefficient 
field does not prevent Liouville’s theorem from failing: Even for smooth uniformly 
elliptic coefficient fields, sublinearly growing harmonic functions are not necessarily 
constant, implying a failure even of the zero-th order Liouville theorem. 

Proposition 21. For any a G (0,1) there exists a smooth, hounded, and uniformly 
elliptic symmetric coefficient field a on R 2 such that the following holds: There 
exists a smooth function u which is a-harmonic and satisfies 

(45) ^ j- u 2 dx'j 2 ~ R a for R 1. 

Proof. By a classical example in dimension d = 2 (cf. e.g. JTD]), for any exponent 
a G (0,1), there exists a uniformly elliptic, symmetric coefficient field a o of a scalar 
equation and a weakly do-harmonic function u o (in particular, it is locally integrable 
and of locally integrable gradient) whose modulus on average grows like |x|“, for 
instance as expressed by 

(46) 

Moreover, in this example 

(47) a o and uq are homogeneous and smooth outside the origin. 
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We now argue that this example may be post-processed to an example of an ev¬ 
erywhere smooth uniformly elliptic symmetric coefficient field a and a smooth a- 
harmonic function u such that still (H51) holds. 

Indeed, because of (l47l) we can easily construct a uniformly elliptic coefficient 
held a that agrees with clq outside of B\ and is smooth. Next we observe that El) 
also implies (using d = 2 and a > 0) that Vu o is locally square integrable, so that 
by Riesz’ representation theorem, there exists a weak solution of 

(48) — V • aS7w = V • (a — ao)Vuo 

in the sense that w and its gradient are locally integrable and that 

(49) J |Vw| 2 dx < (7(A). 

Equation (H51) is made such that u = Uo+w is a weak solution (i. e. locally integrable 
with locally integrable gradient) of 

— V • aVu = 0, 


and thus smooth since a is smooth by classical uniqueness and regularity results. 
It remains to give the argument in favor of (l45l) , which in view of (l46l) follows once 
we show that (14911 implies in particular for large R 


(50) 


(J- w 2 dx'j 2 = o(R a ). 


This is a well-known argument related to “bounded mean oscillation”: By Poincare’s 
estimate with mean value zero we have on every dyadic ball around the origin 

JL 

f (w — -f w) 2 dx j < C(d) ■ 2 n f f |Vw| 2 dx 

J Bon J Bon J \*/ Bon 


which for d = 2 takes on the form 
(51) 


(f ( w-l w) 2 dx\ 2 < C ( f |Vu;| 2 dxY £ C(A). 

Vf / B2 n " B^n J \«/ B^n J 


By Jensen’s and the triangle inequality, this yields in particular ^wdx — 

w dx\ < (7(A) and thus, since we may w. 1. o. g. assume f Bi w dx = 0, 
\f B n wdx\< nC(\). Inserting this back into (15T1) gives 


'/ 

B2 n 


dx I < nC(X), 


that is, (l50l) in the stronger form of 

(J- w 2 dx^j 2 < (7(d) log R. 


□ 
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